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1.0 SUMMARY 

In this study a mathematical tool Is developed for evaluation of 
antenna radiation pattern effects on the Shuttle Orblter S-Band 
Quad antennas. A ray optics approach is used which Includes multiple 
Internal reflections with special consideration to reflection from the 
metallic Orbiter skin. The study shows significant depolarization may 
occur as the angle from the normal increases. Also the effect of 
changing tile thickness on the beamwidth of the upper Quads. versus the 
lower Quads results in a small increase in the lower Quad beamwidth 
compared with the beamwidth in the upper Quads. The results of this 

study may be used to minimize testing in the optimization and evaluation 

% • ■ 

of the S-Band Quads and the computer tool develoDed in this study may 
be used to evaluate other Shuttle Orbiter antennas. 
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2.0 INTRODUCTION 

The purpose of this study is to evaluate theoretically the effects of 
the dielectric materials which cover the S-Band Quad antennas on the 
Shuttle Orbiter. The antenna is assumed to be a point source and the 
ray optics method using multiple internal reflections to calculate 
the tranmission coefficient as a function of angle is used. Two thicknesses 
of material corresponding to the two upper and the two lower Quads 
are considered. In this study a total of 20 individual rays are 
considered to improve the accuracy of that obtained in an interim 
report (1*.2-DN-B0403-004 dated November 21 , 1976). The previous report 
included only 6 individual rays and showed more severe radiation pattern 
effects from the multiple layer dielectric thermal protection system. 

I The results of this study may also be- used to evaluate the radiation 

characteristics of other multiple dielectric covered antennas thus 
minimizing the need for testing to determine optimum design. 

In this paper background information is given describing the specific 
problem and previous work. The theoretical basis for the solution and 
the formulation used are then described. This is followed by the 
computer results and conclusions. 



I 


1.2-DN-B0403-005 
Page 3 


3.0 DISCUSSION 

This section contains (1) background Information, (2) the theoretical 
basis for the ray optics technique and (3) a description of the 
equations used in the computer formulation. 

3.1 Background 

Traditional solutions for dielectric and radome covered antennas have 
involved the use of plane wave transmission theory through multiple 
dielectric layers (References A, B and C). More recent exact solutions 
have invoived the Fourier transform technique for a sing’le dielectric 
layer with an assumed aperture distribution (References D and E), Another 
method considered was the method of moments (Reference F), which 
utilizes mutual coupling to produce the antenna pattern. Both the 
Fourier transform technique and method of moments have been developed 
only for single dielectric cases. Because of simplicity and 
adaptability, it was decided to pursue a modified version of the 
plane wave transmission theory including multiple internal reflections 
and the effect of the ground plane reflection. Previous methods have 
not considered the ground plane reflection. 

The thermal protection system (TPS) consists of Lockheed LI-900 tiles 
having dimensions of 6 inches by 6 inches and varying thicknesses. 

The antenna TPS interface is shown in Figure 1. The multiple dielectric 
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FIGURE 1. MULTIPLE LAYER TPS ANTENNA INTERFACE CONFIGURATION 
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layers are seen to consist of. an RTV-566 layer, Nomex Nylon (felt), a 

second layer of RTV-566, the silica portion of the TPS tile and a 

waterproof coating of borosillcate. The most significant part of 

the dielectric interface is the LI-900 tile which has a thickness of 

1.68-2.15 inches hear the two lower Quad antennas and 0.41-0.12 

Inches near the two upper Quad antennas. The tiles are 

placed on the metal skin of the Orbiter in an alternating fashion . 

similar to a house brick matrix. The spacing between the tiles is .060 

inches and the walls of the tiles are coated with borosillcate. Since 

the thickness of the bcrosilicate coating and the spacing between tiles 

Is small compared with a wavelength at S-Band, the effects of the walls 

are not incorporated in this study. The basic techniques of calculating the 

angles of refraction, transmission coefficients and' reflection coefficients are 

developed in this study and later could be applied to wall effects which 

may be appreciable at Ku-Band and higher frequencies. 

3.2 Theoretical Basis 

This study assumes an isotropic hemispherical radiator from a point 
source with individual rays incident upon five dielectric materials as 
shown in Figure 1. The angles of refraction are determined by Snell's law 
such that 


stne i . _/E. 

sin 0. /Tj - 

J ‘ 


* ( 1 ) 
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where 


Is the input angle of incidence measured from the 
Internal normal in the 1 th dielectric 


e j 

Is the output’ angle measured from the internal 
In the dielectric 

normal 



Is the relative dielectric constant In the 1 th 

dielectric 


c j 

Is the relative dielectric constant in the 

dielectric 


Equation (1) may be rewritten in the form below to calculate successive 

refraction angles such that 

f £ 7 • 

0. = arcsin ^ ■ • sin 9.) (?) 


The tabulation in Figure 2 shows the angles of refraction through each 
layer of dielectric material. The input angle is designated as 0 O . It 
Is noted that the input and output angles are the same since e 0 * e« c 1.0. 
A special case exists when the ray is passing from a material of high 
dielectric constant to one of low dielectric constant such that 
complete internal reflection occurs. This happens when the Incident angle 
Is equal to or greater than the critical angle (© c ) where 


0^ c = arcsin 



(3) 


The critical angles associated with the high dielectric constant materials 


are found to be 



G lc ' 

71.57° 

(RTV-566) 

G 3c ’ 

32.74° 

(RTV-566) 

0 5c “ 

27.15° 

(Coating) 
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As may be noted In Figure 2 npne of the previous critical angles are 
exceeded. If, however, the dielectric constant of the material at 
the antenna surface Is significantly greater than air (e # « 1), then 
complete internal reflection will exist at angles equal to and 
greater than the critical angle. 


The second factor to be considered Is. that of the phase delay which 
occurs In each dielectric material. From geometrical and electrical 
considerations the phase delay in the material Is given by 




2tt f 


/s 


(4) 


COS 0 


n 


where 


f 

c r 

- X r 

C 

G 


Is the phase delay In radians 
Is the frequency in MHz (assume 2200) 


th 


Is the relative dielectric constant of -the n dielectric 

is the thickness of the n th dielectric 

Is the speed of light (assume 11808 megainches/sec.) 

th 


is the angle of refraction in the n dielectric 

li 

The largest phase delay takes place in the silica part of the TPS and 
ranges from 135° to 353.7° for the two lower Quads and from 29.4° to 
77° for the two upper Quads. The incidence angles for the above 
delays range between 0° and 89°. The delay in the fc*lt material ranges 
from 31.8° to 37.4° for the above conditions and the delays in the 
RTV-566 and the borosilicate coating are less than 2.3° degrees. In 
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general, single layer phase shifts of 180°, 360° etc. will provide 
maximum transmission and those of 90°, 270° e c. will provide minimum 
transmission when Internal reflections are considered. 


The next factor considered is the transmission coefficient for 
perpendicular polarization which is the ratio of the transmitted 
electric field to the Incident electric field given by (Reference G) 






cos 0^ 



cos 0. 

J 


(5) . 


where Is the transmission coefficient for perpendicular polarization 

Jj - which Is a complex number 


Is the transmitted field which is perpendicular to the 
plane of propagation . 

A 

^Ij^lnc Is the Incident field which is perpendicular to the plane 
of propagation 

e.| is the relative complex permittivity in the I** 1 material 

Cj is the relative complex permittivity in the j* h material 

0,j and 0- - input and output angles measured from their 
J respective normals. 

The relative complex permittivity may be expressed as 



c • J c. tan 6 
n n n 


( 6 ) 
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where e Is the relative dielectric constant In the n**? material 
n 

tan is the loss tangent In the n th dielectric 


For parallel polarization the transmission coefficient Is given by 
(Reference G) 

A 

T ijn " 


2 /e c1 cos0 1 


(7) 


/F C 1 cos 0j ♦ /T ci cos 0, 


A t fh 

where x^ () is the ratio of the transmitted electric field in j ’ 
medium to the incident electric field in the 1 th medium which is 
parallel. tj the plane of propagation. Other variables are similar 
to those associated with Equations (5) and (6). 


* . th 

The expression for the reflection coefficient r jj x 1n the ^ dielectric 
from the j th dielectric is given by (Reference G) 



Eljj- ref 

A 

^ij_inc 


/*c\ cos 0. - /t ci COS 0j 

cos o, + sr cj cos Sj 


( 8 ) 


where f. . „ is the incident electric field in the 1^ medium with 
nc 

perpendicular polarization and E- Is the reflected electric 
field In the i th med ium with perpendicular polarization. Other 
variables are defined in Equations (5) and (6). 


For parallel polarization the reflection coefficient Is given by 
(Reference G) 
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* . 111 ref . cos e i "'Ti cos 0 i (9) 

1j» ? , 

mine • / e c1 cos Gj +/ e cj . cos 0^ 

It is interesting to note that parallel polarization has associated with 
it an angle where the reflection coefficient becomes zero 
called the Brewster angle or the polarizing angle. 


This angle is given by 



tan 


-1 



( 10 ) 


The Brewster angles for the S-Band Quad dielectric layers are 
found to be 


0 OB 

= 63.43° 

CO 

M 

© 

= 28.41 

© 

CO 

= 43.49° 

' °4B 

= 63.72 

0 2B 

= 46.51° 

0 SB 

* 24.53 


The preceeding refraction angles correspond to reflection angles 
at which the internal reflected field will only be perpendicularly 
polarized for each respective dielectric. 


Another factor to be considered is the attenuation in each dielectric 
material which is expressed for a .low loss dielectric (tan <5«1) as 
(Reference H) 


off /T x n tan (6 n ) 


C cosO in 


(ID 


» 
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where 


°n 

Is In Nepers 
• 

f 

Is the frequency in MHZ 

e n 

is the dielectric constant in the n** 1 dielectric 

X n 

is the thickness in inches of the n th dielectric 

tan (5 n ) 

is the loss tangent of the n^ dielectric 

C 

is the speed of light in megai. hes/sec. 

9 i 

4l| 

is the angle of incidence in the n, dielectric 


The actual attenuation of the electric field is given by 

ATTn = e "n (12) 


It should be pointed out that the power transmission coefficient for 
a circularly polarized incident ray may be obtained from the following 
relationship 


|T|‘ 


L. + T n 

— Z 


(13a) 


or 


|T| Z • 1/4 { |T„] 2 + |T ,| 2 + 2 |T ||T„| cosS > (13b) 


where T |( is the overall transmission coefficient for parallel polarization 

A 

"j_ is the overall transmission coefficient for perpendicular polarization 

is the difference in phase between perpendicular and the parallel 
overall transmission coefficients. 
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3.3 Formulation 

The computer formulation to calculate the transmission coefficient for 
circular polarization (as generated by Orbiter S-Band Quad Elements) 
involves (tensive use of the previously developed formulas. The 
transmission coefficient for a direct ray may be written in the 
following notation for perpendicular polarization 

A A A A A A A * 

T 06 ± = t 01_l t 1«?_l t 23j_ t *34_l t 45_l t 56 ± 

x e" a * e' a2 e~ a 3 e"® 4 e" ay (14) 

x e _ J* (A<t>i + A4>2 + + A4> 4 + A4>g ) 

% • • . 
where A<>, t and a are computed from Equations (4), (5) and (11). 

For parallel polarization the symbol "j_" is replaced by/'n". The 

A _ t 

t's represent complex transmission coefficients, e s represent 
attenuations and s represent phase delays in passing through 

each dielectric. The subscripts designate the following 

0 for the antenna cavity 

1 for the inner RTV bond 

2 for the felt (strain isolation pad) 

3 for the outer RTV bond 

4 for the silica portion of the tile 

5 for the borosilicate (waterproof coating) 

6 for the air or vacuum region outside the borosilicate coating 


I . I I I l l t 
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Since there exists the possibility of an Infinite number of Internal 

reflections and subsequent retransmissions, the transmission and 

reflection coefficient magnitudes for each interface were determined 

In order to evaluate only the most significant reflections and retransmissions. 

A matrix showing some of the data is given in Figure 3 for a total tile 

thickness of 1.87 inches with the weatherproof coating. The 

reflection and transmission coefficient magnitudes are included as 

well as the phase delay in degrees and the electric field attenuation 

factors. It is observed that some transmission coefficient 

magnitudes exceed 1.0 in going from a material of high dielectric 
• • • 

constant (low impedance) to one of lower dielectric constant (higher 
Impedance). An analogous result occurs when a transmission line of 
low impedance is connected to one of high impedance in which the 
transmitted voltage may double in magnitude in the high impedance line. 

From Figure 3 it is observed that the reflection coefficient 
magnitudes for r^ and ^3 are quite small (less than .04) and that 

A A A 

the reflection coefficient magnitudes for r 34 , r 45 and are 

somewhat more significant being on the order of .3 for 0 = 1°. Internal 

# . /N /\ 

reflections resulting from r^, r^ 5 and r^ are shown in equations (15), 

(16) and (17). 
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FIGURE 3. TRANSMISSION AND REFLECTION COEFFICIENT MAGNITUDES AND OTHER DATA 
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R56 1 


R45, 


T 01_L T |2_l T 23 jl J34jl T 45j_ T 56 ± 

A A A A A A 

x t 54j_ t 43j_ t 32j_ t 21j_ r IOj_ T 12j_ 


05 ) 


A A 


-3 (a ] + a 2 + a 3 + a 4 + 


x T 23 i T 34 L T 45 1 T 56 ii 
x e -J3(A<J> 1 + A4> 2 + A<t>2 + A$ 4 + A<f> 5 j 

A A A A . A A AAA 

t 01jl X I2j_ 17 23_l ^ 34_l r 45 jl T 43 J _‘ t 32_^ 21 ± r l 0 


A A 


AAA 


X X T 12_l T 23 ± T 34j_ t 45j_ t 56jl 


-3( a + a + a + a ) - a_ 

xe v l 2 3 4' e 5 


x e 


•j 3 (A$.| + A<}> 2 + A4> 3 + A^ 4 ) e -jA(J) 5 


06 ) 


R34 


T 01j_ T 1 2j_ t 23j_ r 34_L T 32_ L T 21 ± r 10j 


X T 01 x 1 42j_ t 23_j_ t 34j_ t 45j_ t 56 ± 

xe^l +a -2 +a 3)e-N; a 5 ) 

j3 (A^ 1 + A4>2 + A^ e' j X + A< ^ 


(17) 


x e 


where r, n = -1 for reflection from the metal skin of the Orbiter. 

IUjl 

Similar expressions for T R2 3 and T R ^ 2 may be written and the above 
equations for parallel polarization may be obtained by changing the 
subscript "J." to "u". In addition to the preceeding internal 
reflections additional reflections are considered because of the 
significant magnitude of the reflection coefficients in layers 3, 4 and 5. 
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These include two rays which reflect from the boros ilicate/air (or 
vacuum) interface r^g and from the sil ica/borosil icate interface 

A A 

r^g and which are then reflected from the si 1 i ca/RTV interface r^ 
and transmitted forward. The equations for these two rays are 


A A A A A A A 

T 546 a * T 01_ l T 12j_ T 23j_ t 34j_ t 45j_ r 56j_ 

(18) 

x t 54j_ t 43j_ t 45j_ *56j_ e ‘ * a l + °2 + °3* 
x e - 3 ( a 4 4-o 5 ) e -j (A*, + A * 2 + A$ 3 ) 

x (A* 4 + A4 5 ) 


A A 


AAA 


T446 




T 01j_ T 1 2j_ T 23j_ T 34 r 45 r -43 


A A 


-(a. + a, + a, ) -3 b, 


x t 45 t 56 e"'”l ’ “2 ’ "3 ' e 4 


x e 


• a 5 e -j (44, + 4« 2 + A* 3 ) e - j 3 aa 4 


(19) 


:■ e' j ^5 


It is noted that r 


43 


= - r 


34 


and that the expression for parallel 


polarization may obtained replacing the "iJ 1 subscripts with "II". 

Similar expressions may be written for ^ 555 ^ which represents the transmission 
coefficient for a single internal reflection in the 5th layer as well as 

^255_l anc * T 355j. W ^ 1C ^ re P resents two and three internal reflections 
in the 5th layer. Also, expressions may be obtained for secondary 


! 
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reflections of three largest reflected rays T R56 ^ T^ and 

A A A 

Tr 34 ^ from the 3/4, 4/5 and 5/6 Interface as T^g^ ^X45j^ anc * 

Tx 34 The last three transmission coefficients represent a total 
of 9 rays. 

The composite transmission coefficient for 20 rays with perpendicular 
polarization becomes 


T 06j. + T R56_l + T R45j_ + T R34l + T R23J_ + T R12 ± 

A A A A A 

+ T 546_l + T 446j_ + T 556_l + T 255i + T 355_l 


+ T X 56j_ + T X4 5 _l + T X34x 


Parallel polarization is obtained by replacing the "J_" symbols with 
"H". The transmission coefficient for circular polarization 
becomes. 


Tx + T„ 
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4.0 RESULTS 

This section describes pertinent results from the computer runs. A 

plot of the direct transmission coefficient magnitude for both 

perpendicular and parallel polarization Is shown in Figure 4. These 

factors do not consider internal reflections and are presented only 

as an aid in understanding the complete effects of the multiple layer 

dielectric covering. It is observed that the perpendicular polarization 

coefficient is significantly smaller than the parallel coefficient 

which shows a tendency for the outgoing wave to have a predominant 
« • • 

parallel polarization. If the input ray is assumed to have 


perfect circular polarization the axial ratio of an outgoing ray may 
be computed by the formula below assuming 6p - 0. 


A. R. 


20 log 



( 22 ) 


At an angle of 50° the hypothetical axial ratio .would be 3.77 dB. At 
70° the hypothetical axial ratio becomes 8.7 dB . 


The results of the computer run with the complete transmission coefficient 
including a dB factor for circular polarization are shown In Figure 5 and a 

plot of the results is given in Figure 6 for both the upper and lower 
Quad antennas. The dB factors represent pattern changes which would 
take place but do not necessarily represent a loss of energy since the 
pattern is redistributed and since energy reflected back into the antenna 
aperture -j s no t specifically treated. 


TOTAL DIRECT TRANSMISSION COEFFICIENT MAGNITUDE 


i 



FIGURE 4 
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FIGURE 5. COMPUTER RESULTS FOR COMPLETE TRANSMISSION COEFFICIENT MAGNITUDES 
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It Is observed from Figures 5 and 6 that the TPS has a noticeable effect 

on the lower Quad antennas with a slight reduction in gain on axis (normal 

to the antenna surface) and relative increases in gain off-axis with a 

maximum at 0 = 45° and 0 ■ 78°. The general effect of the TPS on the 

lower Quads will be to slightly increase the antenna beamwidth over that 

obtained when the antenna is operated without a TPS covering. The effect 

of the TPS on the upper Quads is to gradually decrease gain as the angle 

off-axis increases. From 0 = 0° to 0 = 45° there is little change in the 

level; however beyond 45° the gain level falls off significantly. The TPS 

in this case has the effect of narrowing the beamwidth of the upper 

Quads. In addition to beamwidth variation effects, the transmission co- 
• * • 

efficients for parallel polarization in Figure 5 are found to be generally 
larger than those for perpendicular polarization. The result is that 
axial ratio degradation occurs off-axis. Using Equation (22) the 
approximate axial ratio degradation for the upper Quads is found to be 
2.9 dB at C = 50° and 6.36 dB at 0= 70°. For the lower Quads the axial 

ratio degradation is found to be 2.25 dB at 0 = 50° and 7.77 dB at G = 70°. 

If the S-Band Quads have a lower parallel polarization components with no 
TPS the previous axial ratio numbers may improve. Also, since t.ie 
maximum and minimum for perp ^dicular and parallel polarization do not 
occur at the same angles some unusual axial ratio variation may be 
exp?ctr J . 


I 
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5.0 CONCLUSIONS 

The developed computer program shows that slight changes in beamwidth 
may be expected between the upper and lower Quad antennas caused by 
the difference in thickness of the TPS tiles. Also, a general degradation 
In axial ratio is shown as the off-axis angle of the Quad antenna is 
increased. The axial ratio degradation may be improved by using an 
antenna design with a less pronounced parallel polarization component when 
operating without a TPS cover. 

This investigation shows that the multiple layer dielectric covering 
for the upper Quad antennas will slightly decrease the antenna beamwidth 
of that obtained without a TPS covering. In addition, the lower Quad 
patterns are slightly flattened on axis and the beamwidth is slightly 
Increased over the no TPS condition. Since the Dattern modification 
effects of the TPS are sensitive to changes in dielectric constant, 
material thickness and frequency, a change in any of these parameters 
will result in a different pattern modification effect which may be used 
to optimize the antenna coverage. 


I 
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LISTING FOR PROGRAM WHICH 
PRODUCES DATA IN TABLE 
FORM (SEE FIGURE 3) 
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K?3pEns(K2*C0S(T2r)-R3*C0S(T3R| )/(R2*C0S(T2R>*R3«C0S(T3R) ) 

— k-jMPEns |kj*lOsIIOR>-RS*COSITMR))/(R3*COS«T3R)^RR»COS(TRR)) — 

KmSPEn=|KR*C0S(T‘I.R)-rs»C0S(TSR))/(RS*CO5(ThR)4-k5»C0S(TSPI) 

M,6PEN»(R-i«COSlTi>f<)-R6*COsUZ»R))/IR»*COSlTSR)*R4*COSCT6Rl) 

C R i 2 p A R - R S 6 P A R Rc.PRc.SlnT COMplLX ELECTRIC FIELD REFLECTION 

~c ~ coefficients for parallel polarization 

Kl2PARc(Ki*COb(T2K)-NZ*COS(T!R))/lRl*COS<T2R)*R?*COStTlR)> 
K?3pAi<=(K2*C0S(T3R4-r?3*cOSIT.'»R| )/(K2«C0S<T3R )-»R3»COSIT2R) ) 
K3iF’ ARa(K3*CUSITMi<)-RH*COb(T3R) )/(KJ*CUSlTMR ) ♦ 9 *» • C 0 S ( T 3 R ) ) 

* Ri4SPAK*(KS»COS(TbR)-RS«CUS(TRR))/|K4«COSnSR)+R5»COS(THR)) 

Nm.pAK3|KS*C0S(T6R»-Ro*C0SlTSR))/|R l j«C0SlT6R)+R6*C0S(TS p )) 

C TSSPLn - I 2 1 P L N HcPKt5c.Nl TRANSMISSION COEFFICIENTS FOR REFLECTED 

c Rays ■=* i r h perpendicular polarization 

I c,NPEn*2» • K S • C 0 $ I T SR I / ( RS • COS I TS/x ) ♦RM«COi> l T M R ) > 
1*43pEn»2**km«COSI IHR)/(K4«CUS(TSK)*R3*COS(73R) ) 


OK1GUNAU tAGBja 

^a^ftQD ATT1T JTlfl 


/ 
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I 

* » 


-C“Ts*iPao 

c 


T32PtN»2»*R3*C0SU3RI/<K3*L0S(T3R>«’R2*C0S<T2R)) 
T 2 lPCN* 2 ««K 2 *C 0 Sir 2 K)/IR 2 *C 0 SU 2 HMN|«CUS>(TIR)| 

AO --I 21 PAK KLPKtstNT TRANSMISSION COEFFICIENTS FOR REFLECTED 
KAYS ..ITH PARALLEL POLARIZATION 
Te,HPAK»2««KS*CUS<rbK)/(NS*LOS(TSR>*'R‘«*COS(T5R)| 
TMipAK*2**K4*CuSUHR)/(RM.LUS(TJRMK3«C0b(TMR)) 

”■ I12PAr* 2« *K3*tOS » T3K I / ( K3*C0S ( T 2 H I ♦R2*COS I T3R ) ) ' 

r?IpAK*2.*K2«C0S< r2R)/IK2*CUSU|RMRI»C0SlT2R» ) 
TrilPt*LAbb(rClPEN) 

ri2pE*LAbSlTl2PE f O 

— 1 ?JPt.«CAOb t T2 JPe.N| - . — 

T3HPt = CAb3(T34.‘ , E'U 
IhSPcsc AB b l T HSPEN ) 

It,oPt = LAbaiTS6Pt , '*» 

— FntPAsLAbb ( T21PA* I 

T | 2PA iC Abb t T 1 2 P A« » 

T ?JpA = L Abb ( T 2 3P A R ) 

T 3 S P A s C A b o ( T 3 S P A K > 

— \ cj 5PA = L Ab5 t T HiP AA 1 — -- - 

Tc,6pA*(. Aba ( T b 6P AR ) 

K | 2PE sC AbS (rt 1 2PE-U 
K?3PE*LAOj l R 2 J P E'< ) 

- ~H 3 R P E = c A & i ( K 3 1 p E <t I 


RHbPE*CAbbtRH 3 PtN) 

KfjfePtsQAbi l R 5 fc P E N I 
K|2PA=LAbjlRl2PAR) 
n?JpAswAb3tf(2jHAi'l 

-K-jHPA^L Abi IR34PAR ) * " 

K m 5 P A s C A b a l R R S P A * ) 

KSoPAscAbalRboPA.w 

rs*iPE = CAbS l TbSPtf* ) 

f tr3pEsCAbb U R JPEN f 

T i2PEsLAbb l T32PEN ) 
r2lpE*LAb.»(T2lPENi 
KSpAsLAbs H&RPAR ) 

f nbpA = LAb^tTR3PAR I — 

[ i2PA*LAbS(T 3 2 P A •< ) 

T? I P AsC Abb l TZI PAR ) 

' AFTT s (At!*AL2*XEi*AEM*X£S)/KD 

at fsCMPLA i C.O .-Xt T T > 1 — ‘ " - 

TpLN=lJlPcN*T12Pi-N*T23PEN*r3HPEN*TH5 l , EN*T56PEW#ATTl*ATT2»ATT3 

l*ATT5*CLXp»xEr) 

Tpt=CAb5(TPEN) 

r p AR * IJiPAR*112PAf%»T23PAR*l3RPAR*TRbPAR*T5toPAR*ATTl*ATT2»ATT3 

1 • A T T 5 • v. E X p i X t r ) 

IpAsChdSI TPAR ) 

«.<ITE(e, l Ib)TCIPEiri2PtiT23PE.T3HPi:.TMsPEiTb&PE 

«0 ITEle> l lt,ML’lPAiTl2PA,T23PA l T3HPA»TRbPAiTS&PA ~ 

nylTLl6»lb>Rl2PE»R23PE»R3sPE»R t tbPt-»P5bPE 
i< N 1 T E lb i 1 a)!< 1 2 P A ir\23PA f K J R P A t R R 3 P A i R b 6 P A 

rtpnElb,lV)T2lPL>;32Pt,TH3PE,TSMPE 

n* ITtlb.»V)T2lPXiT32PA,TH3PA,r5SPA : 

iniITE(6,ld)ATTli*T12|ATT3 l ATTH t ATTi 

u o l T E l 6 , 2 3 n P L 
n * 1 T E ( o , 2 j ) T P A 

lb Fn«ttAl ( lOx ,FS.?,b USX ,Fb. 2 ») • ' "* ‘ " 

lo t unMAT l 3 ua if- 3*2 1 1 t 1 SX »F 2 )) 

19 FoRMAl (30x.F3.2,3l 15X.FS.2J I 

2b FoKMAUl2j>A,Fb.H) 

20 £ u in T I i« J E 

J E J ♦ I 

IrU'Lti.l )■ GO T0‘ 3 ■ ~ ----- 

lnJ 
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LISTING FOR PROGRAM 
WHICH INCLUDES MULTIPLE 
INTERNAL REFLECTIONS 
TO PRODUCE COMPLETE TRANSMISSION 
COEFFICIENTS (SEE FIGURE 6) 


I 

I 


1 


REAL IE 1 

— COMPLEX 

COMPLEX 0 0,01 ,R2 ,PJ 


l ,IE2, IE: ,!E“ , IE5 
T P E N , * E T , T o t q 


complex 
complex 
complex 
COMPL E x 
complex 

COMPLEX 

COMPLEX 

COMPLEX 

complex 
complex 


COMPLEX o?.,oi,0?,cj.3P,0'.R6 

COMPLEX TC1 p E*4,t 1 ? 3 f N , r?2PEN,T3<*PE**,T«5PEN. 
COMPLEX TE.lPio.Tiip/o.TZJPAR.TjMPaP.TUSPAP, 

COMPLEX Oi:PEN,O^TPr f , t ;j(.PCN,345PEN,P p bO£N 

COMPLEX RlZPi= , ,R:;P.“?.»34PAP,P4;PfiS.',PE e PAR 
COMPLEX EC3 .EC!,EC?,EC:,ECh,ICS,EC 6 
T54P*P t T“’P*P,T32PAP,T21PAP 
TSUPC*:, Tu'sfn f T32P£N, T21PCN 
TPE »,0 T, TPAS f T , TC 1® 
TP5EP r ,TPCfc p A,TOn<;PE,Tf>««5PA,TR3‘4PE, 
T54<;PE f T;-tPfl,TUi<6P£,T4»*6 p A 
TE56PE, T bii' i 


X t 4fc , X44 o,X55o 
xetp4e,xet?:l 
XCl.XCZ 


c 

c 


W vn r r a 

COMPLEX XC1 ,XCZ, XCJ .XC4.XC5 

L-COMPLEX T 2 R ? P E ,T ZEE Pi , I 3 5 N P £ , T 3 5 5 r A 

COMPLEX TR?:PE t TP?3PA,TPlZPE,TRirPA 
COMPLEX TX34PE,Tx34rA,TX45PE,TX4SPA,IX56PE, 
J - C 

THICKNESS E OR UPPER CUACS - 

xu:.405 

GO To 4 

THICKNESS EOP-LOwER LOADS 

3 X4 = 1.S6 . . 

* WO I T E (6,5) 

5 FORMA niHl,lCX,5HTNfTA,6X,3HT'’E,7*,3HTPA,6X 

XI - X5 REPRESENT DIELECTRIC LATER THICKNESSES 
X12.015 
XZ=.2 C 
X3U015 

- xsuoi - •. 

RI = J. 14155 

EO - E 6 REPRESENT RELATIVE DIELECTRIC CONSTANTS 
EC=1. . 

r « - ii m 


C TDl 


El = 4 

# "1 

E?- 3 

.6 

E3 = 4 

# n 

E«»U 

• 17 

E5-4 

, Q 

E6U 

.0 

- TDS re 

TDl : 

.CDS 

TD2 = 

.04 

TD3: 

. CO 5 
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T56PEN 

T56PAR 


TR 34 P A 


T X 5 6 P A 


,*»htc io, 6 x , eh tcir ( re ) / ) 
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T0ls.r016 
TP5 = . CC 3 

-C—lCl - I E 5 REPRESENT IMAGINARY PARTS-OF-RELATIVE COMPLEX PERMITTIVITIES 
IF 1 :-E 1 • T? 1 
IE2=-E2*TC2 
IE3=-E J*TC3 

IC*«z— r«*«TD** . .. 

IE5 = -E5» T25 

C ECO - E C fa REPRESENT RFLMIVE COMPLEX PERMITTIVITIES 
ECOiC-PLx ) 

ECUCl'PLXIE Itl-l ) * — — .... 

EC2 = C“PLX I E 2 , ! E 2 > 

EC3=C“PLX IE T, IE ! I 
EC«-C“PL* ( E u , IE4 I 

ECS = C “PL x ( E c , I E 5 > - . . . 

EC6=C“PLX (Efa, r.O ) 

C RO - Rfa REPRESENTS THE FELATIVE REFRACTION INDICES 
RP=CSCRT JEC") 

RISCSORTIECH 

R2 = CSC»T (ECD 
R3=CSORT (EC 3 I 
RM = CSCRT ( E C *•' I 

RS-CSCRT l 2 C 5 I . _ 

Rfa = CSCRT (ECfa ) 

C C IS THE SPEED OF LIGHT IN megain-ches/sec 
C= lloD8 . 

-C F IS FRECUEfJCY I.N MEGAHERTZ _ 

p ■ 2 ^ ^ ^ 

C RO CONVERTS RADIANS TO CEGREES 
RC=16D./3. 1M5<? 

DO 20 1 = 1,35 . • _. 

TC = FLOAT ( I » 

C 

C TOR - T fa R REPRESENT ANGLES OF REFP ACTION IN RADIANS 

; T3/RD 

TlRrASINfsORT tEC/F 1 « TC° ) » 

T?P=a SIN I S : ~ T (: 1 /E2 ) • $ IN < T 1= ) > 

TJR = ASIM S.-RT lE2/ r !)*SI* l T^R ) ) 

TURZASIMSC^T IZI/EV J#SIN(TI- ) ) .... . ... 

T5R = ASIMSCRT(EV/ES)«SIMThR)I 
TfcRrASIMSOPT(E5/ES)*SIN(T5=»l 
C TO -T6 REPRESENT ANGLES Or REFRACTION IN DEGREES 

- — . . T1 = T1R*RC ........ 

T2=T2R»R0 
T3=T 3 p *PC 
T<i = T*4R«RD 

TS=T5R*RD .... .. . 

Tfa=TfaR«R: 

C X E l - X E S REPRESENT PMA*E C r L a y S IN DEGREES IN EACH OIELECTPIC 
XE1=S0RT IE 1 I • 3faC*F*Vl / (E»CC C ( T 1 P I ) 

XE7 = SCRT IE2J*JiC»F»x«;/lC*CCSl T 2 R 1 ) .... 

XE 3 = S OR T (f 31* “ *F • * I / ( C • C c S I T2R | ) 

XE**:SCRT (• « l • >* «r *»m / if.c CS( T u p ) j 
XE 5 = SCRT • E 5 1* 3fc r vF • X S/ ( C • CCS ( TER ) ) 

C 

C XRI - XRS REPRESENT PHASE DELAYS IN RADIANS IN EACH OIELECTPIC 
XR1=XE1/RC 


OKI 1 

ov 


POOR Q^aVUTY 



1 


r 


;] 
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xR2=xE2/ p n 

XP3-XE3/ p T 

XPn:XEu/£>D - - 

XP5:XES/RP 

C XC1 “ X C 5 REPRESENT PHASE DELAYS IN RADIANS EXPRESSED AS A CNPLX NO. 
XC1 = C “PL x CC.u ,-*»•! I 

XC2 = ChPlx I 

XC3=C“ p LX (3.r,-XP3) 

XC« = CfPLX 10. j.-xRM 
XCS:C“PLX ( j.C ,- xp= ) 

-C- ATT 1 - A T T j REPRESENT ELECTRIC PIELO ATTENUATION FACTORS 
ATTlrEx p (*l • « (PI « r • ^ C * T ( E 1 )«T A N ( T T) 1 ) * x 1 ) / (C*COS( T 19) 1 
AT72:Ex p f-l.*(RI*P*SCPT(t::)-TAMTrZ)«xr)/(C*COS(T-'PM 
ATT3:ExP(-l.*(PI»r*SQ>T(£3l«TfN(Tr3l»x3)/(C»COS(T:S)) 

ATT<*rExP(-l.“(PI«‘f*SC-T{E“)*TA.N(T''‘*)»xu|/(C*CCSCT«P|j 

ATT5:ExPJ-l.*(PI« r *'C°TtE=)*TAN(T r '5)cX5)/(C»COS«TSP)) 


C TC 1 PEN - T 5 fcP E N REPRESENT CCf' p LFx ELECTRIC P I EL P 
C CCEPFICIENTS POP °ERPENO ICUL AR PO 

. TC- 1 PENR 2 . •RP»COS I T 3 P ) / ( p C«COS ( TOR > *R 1 -COS I T IB 

T 12 PEN: 2 .*Rl»C 0 S(TIC)/(P:»CCS(Tl?)«cr«C 0 StTP» 
T 23 PENR 2 . * P ? • C 0 S l T 2 F ) / IRE-COS ( T 2 R ) *P J • CO < T ’ R 
T 34 PEN = 2 .*R**CCSlT 3 F>/(R 3 *LOS»T 3 W)*R 4 -COS(T«<P 

TN 5 p ESR 2 «* pt i*C 0 S ( TUP > / | R •« « L C S I TuSJ ♦ R E * C 0 S I TSB 

T 56 PE*.: 2 .»PE*»C 0 SlT?PJ/(PE«C 0 S(T 5 B)*=s-CCi(TO 
C T 01 PAR - T 56 PAB REPRESENT CC“ p LEX ELECTRIC FIELD 
C CCEPFICIENTS F C p PARALLEL PpLAPIZ 

T 01 PAP: 2 .*R n *C 0 S ITCRWlRC-CCSlTlRJ.Rl.COslT'JJ 

T12 p AR:2.«Rl*CCS l TIP ) /IP 1 • C C s ( T 2 B ) ♦ R 2 “ C C S I T l R 
T23 p A p :2.‘R?»COS(T:P)/(R t «CCS(T3R»*RI’COS(T'P 
T 3 f *PAPr 2 .*Rx«'-CSlTJP)/(P:«.c:S(T<iP)»R‘*»C 0 S(T 7 R 
TuSPARr 2 .*KU»cCSJTuC)/(?usCC 3 <T 5 B)»RS*L 33 lTup 

TS 6 PAPr 2 .»REaCCS l IE® ) / I R 5 0 C C S ( Tt>R l ♦ 3 t * C C E I TER 
C P 1 2 p E N - R 56 PEN REPRESENT :C“ p lEx ELECTRIC " I r L T 1 


coefficients ®eppeno!Clla? 

Rl?PEN=«P!*C0SIT!R)-B2tCCS<T2R)»/(R:»CCS«T!R» 

R23PEVS(R?*CCSfT23|-$**ro«|T3 = ll/<R“*fCSlT? e l 
P3APEN:(R?«CCS(TIP)-P4wCCStT4F»)/(Rl*CCSI T 3=) 
R4 5PEN:«R‘'»COS(T,P)-RE»r;riTS p n/«P‘*eC:3iT'MP> 
R56PENA(R;*C0SlTSRl-Rfe«CCS»T6RJ)/(PE*C0S(TS p l 
R12PAP - RS6PAR R rp? cSE‘ T CO* ,p LEX ELECTRIC FIEL r ' 


COEFFICIENTS PARALLEL ROla-'TZ 

RirPARrCR^^CSlTEPI-crtCCSt'. l^JJ/lFl-CCSlTERl 

R22PAP:(r?r*CrSlTiP)-K;»C0NlT25U/!K'0 •CC f <TJR> 

R34PaRa(RI»'*'SIT-P)-Ra»'-(.‘. tTIRII/IAltCJSI^AI 
R«SPA p :(R4«rcS(TECj-PE»C:S(Turjj/(t4*<:cSl T;? ) 
RS6FA p :j pt »CCS(T 3 P)-R£-v-CC^ ( T S *? I )/l*?c«ccS ITS®) 
-C-T54PEN - T2IPEN REPRESENT T - - ■ , - I S S I C N CCEFFICIEN 
c Cite WITH r-.A u N T’.'LiF pcla- ::at 

T54PENR2.* pe *C0S (TE = I / ( F - « . ' S ( T S A > ♦ P « « C C S I Tub 
TN3PrN:r.»R‘-.::scTuP>/(Ct.«cos(T«iRi*s:tco'»«TTp 

T 32PE N;2 . *R 3*CCStT2F>/(R:«C0i{I3n)*R2*C0SCT’R 

T2iPFMr.« ; ?:»cc:tT2r)/t !: :“cos(T:R)*Ri*cosiTit’ 
C T54PAP - T21RAR PEPRlSENT T A N S “ I S S I C N C. r FF;CI r N 
c f < v ' with parallel p r lap !eat;:n 

I54PAP = r. *ccs lT' c )/(RL->cCSlT4P|«Ri.*CCSlTSB 

T43PAP;2.*R4»CCS ITUPI/I^-COS ( T 3 R l-ftl-COilTui? 
T32PAQ=2. •R2-CCS C T3P >/ (PJ-COS ITER! -P2*COS « T IP 


RAN SCISSION 
ARI2ATI0N 

RANSHISSION 

TICN 

' \ 

I k 

> \ 


E p L r CTION 
POLARIZATION 

R2-C3S ( T TR l ) 

R I • r C « ( T 3 o » ) 
sj **CCS(T4p;i 
°5 • 0 0 S < T 5® i) 
°F*COS (T 0 R II 
EFLECTION 
TIC*. 

»?-COS l T JO » 1 

P T »COS( TEC i ) 

B4 .COS 1 T 3® )) 
p *j*CC< I T-R ) # 

®6 ** CO S I T S3 ) ) 

S FC® reflected 

CN 

I 


FCR REFLECTED 


nonnn noon . nnoh ono r»o non noon | nno ooo 


! 
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121 PA Fr 2. *R2*C0S IT TP l/(R2*C0S (UR I ♦« 1*C0S«T ?R» » 

C 

-C--XETT REPRESENTS PHASE DELAYS If. RADIANS THRU 5 TPS MATERIALS 
XE7T=lxEl»xE:«xEi*xf4«x£5l/RD 
XET=Chplx (D.C.-xLTT ) 


T PEN REPRESENTS DIRECT TRANSMISSION COEFFICIENT FOR PERPENDICULAR 
POLARIZATION 

TPEN:TDlPEN«Tl2PEN*T2!PEN*T3‘*PEN*Tii5PEN*T56 t> EN*ATTl*ATT2*ATT3*ATT4 
1 *AT TS*CE /P ( XE T ) 


TPAR PEPR r SENTS DIRECT TRANSMISSION COEFFICIENT FOR PARALLEL 
POLARIZATION 

TPAP:TaiPA»*T12PAR*T22PAR*T34PAR*T**5P»R*T56PAR*ATTl*ATT2*ATT3«ATT* 
1 «A 1 1 S *CE xp (it E I) 

TR56PE REPRESENTS TRANSMISSION COEFFICIENT TOR RAT PEFLECTEO FROM 
THE 5/6 INTERFACE AND FRO" ThE hETAL GROUNO PLANE wITh 
, perpendicular polarization 

TR56P r =TClP r ‘.»Ti;P r, *Tr3P c 'N»T3<4PEv»T4SPEN*ost)PrN»T5«iPrN*TuJPEN»TJ? 
1PEN*T2IPEN'»-!.*T12 S EN*T22 P EN*T }4 f EN • T 4 5P E N* T S6PE N* ( (ATT1*ATT£*ATTJ 
2*ATT4*ATT5)«»(2.-> ).CE>P*2.-xET» . 

_I P. 5 6 F A REPRESENTS TRANSMISSION COEFFICIENT PCS °A Y REFLECTED FROM 
THE 5/6 INTERFACE AND FRO" ThE f*ET«L GROUND P L A NE WITH 

parallel polarization 

TR56PA:TCIPAR»TlZ c AC«T::PA-»T34PiP*TtepAR«c>5fcPAR*T t ;upip*T43PAP*TI‘* 

lP4R*T2lPA' , *-l.*Ti:?iR*T2IPA3«TSMPAR*TMSPAR*T56PAR*(iATTl*iTT2»ATT3 

2 *AT T4 • AT T F ) *«» ( 3 . J ) I • CE X° ( I . •* E T ) 

XTR45 = (3.*lxEl*x ! .r»*E;*xE , *l*X':5)/PD 
XLTR4S REPRESENTS P U «SE CEwAy In RADIANS c C c REFLECTION F 0 C *’ .ThE **/S 
. . INTERFACE AND R-FEFLECTICN FRCf THE Ml £ ’ A L ORCuNC PLANE 

XETR4 c ;Cmplx (C.O , - » T9u>) 

TR45FE PE p PES r NT« T - an S " I Si I CN COEFFICIENT F’CC »it REFLECTED FRO* 

The 4/5 INTERFACE a *.C ThE m.ETAl GROUND PLANE wITh 

PERPENDICULAR POLARIZATION 

TP4 5PF : T D pc N • T : . I? , rr...::u 5 .-Ii« 3 ucpr » 1 «Tujt>r», « r t;pr >,* jp ] ?f >j». i , 

l4T12PFN*Tr3PEN»TlaRrN«T45REN*T5bPEN»((AfT!*AlT?*ATT**ATTui«*(2.n* 
2ATT1*ATT2 «aTT 3*t TT4 <.aTT5»CExP ( a t T - ** 5 ) 


TP45PJ REPRESENTS TRANSMISSION COEFFICIENT F 
ThE <* / 5 IMEFFffF AN' R = 0" Ihi. METAL 

parallel R : i a - ■ . • 7 ; ; 

TR45PAITC1P*;.: : . ■ tCaTZ?PAh«T3«PAR*S45PA 

l^TliPAR •TZJf'AR*’ ’i.?fC »Tu;r A3* ^cr’fA* I (AT 
2*ATTl»ATT:»ATT2*ATTi.eMTi*C r xR(xETau5l 

XTR34r(3.«(xEl*XL2»XE'»l«*Ft.*XE5>/R0 


C p 3|* pe c LEC t ED FROM 
oROUND D LANr mITh 

R *Tu»P1P»T'»7PAR*TZ1 p A c *-1. 

1 i»AIT2*ATT».«iT4»»*( ?. ) I • 


XE TR 34 REPRESENTS PHASE DE L * Y IN RADIANS FOR °E r LEC T I ON F D Of THE 3/4 
INTERFACE AND P E r E 1 F C T I ON FPC" THE m£TAL GROUND PLANE 
.... XETR34;C" P LX ( r.C ,-x TR 3s) 


T R 34 PE REPRESENTS TRANSMISSION COEFFICIENT FC® RAY P E p L E C T E T F R C H 
.. _ ThF 3/4 INTERFACE AND FRO." Iht N.ETAL GROUND PLANE WITH 

PERPENDICULAR r c L A ? I Z A T 1 C *. 

TR34PE-TCl t, EN»T12P[N*T23F’EN*RJ4PEN»T3 2PEN«T?IPEN*-l.-*T12PEN»TZ!F’£N 
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1*T34PEN*T45PEN*TS5PENy<»ATT1*ATT2*ATT3>«*«2.»>*ATT1*ATT2*ATT3*ATT4 

2MTT5*CExPlUTfmi 

C 

C TR34PA REPRESENTS TP A NSm I S S I ON COEFFICIENT FOP pay REFLECTED FROM 
C T*f 3 / u I N T F ® F ACE A\0 FCCM THE. mETAL GROUNO PLANE wITm 

C -■ PARALLEL POLARIZATION 

TR3t4Pe:T'!jP;P*Ti:P« c » T ?!PAL.R3aPAP*TT?p 4 R,T?lPAR«-l.*Ti;pAP*T?3PAP 
l*T3‘tPAR*Tu5 :> :P'‘TSt.P«R*(iATTl*ATT2*ATT3)**<2.>»*ATTl*ATT2»ATT3*ATT«» 
2«ATT5*CEXP( xETRJn > 

C 1R23PE REPRESFN T S TRANSMISSION COEFFICIENT FOR 0 A Y REFLECTED FROM 
C The 2/3 INTERFACE a\D From THE METal GROUNO PLANE WITH 

_C - PERPE NCI CUL AR POLARIZATION 

T023P c :TPES*T12PLN*P23PEN*T21PEN*«-I.I*1IATT1*ATT2)**2.)*CEXPI2.*I 
1 XC 1 ♦ X C 2 ) » 

c 

C ~TR23 PA RE PRE SENT s~ "tr”aNS *'1 sYl CN~ COEFF I Cl E N1 FOR PAY REFLECTED FROM 
C ‘ THE 2/3 INTERFACE AND FROM THE “ETAL GROUNO PLANE WITH 

C PARALLEL POLARIZATION 

TR23P A i TP AS «I 12P AR*R23?AR*T2IPEN*I-1 • A « C CAT U*ATT2>**2. » *CEXP 12. •! 

1 XC1*XC2 I ) 

C 

-C TR 12 PE REPRESENTS TRANSMISSION COEFFICIENT F C R PAY REFLECTED FROM 

C THE I/O INTERFACE AND from THE "ETAL GROUNO PLANE «ITh 

c plppen:-icwLacp:lapi?*tidn 

TR12PE:TPEN*p12PEN*«-I.J*ATT1*ATT1*CEXP(2.*XC1I 

f C TR12PA PE°RESE'«TS TRANSMISSION COEFFICIENT fop pay PE r t E C T E 0 FROM 

C ThE 1/2 INTERFACE J ‘ 1 FROM The METAL GROUNO PLANE -ITH 

-X PARALLEL PCL A A T 2 A T I <,« 

1P12PA=TRar*R120as*(-i.I«aTT:«aTT!*CEXPI 2.*XCII * 
XW5Hfa:(XEl*xEr**D3*:.-»*E‘-*xE5n/ 0 t) 

A546rCMPLX«n.a,-x-E4 D > 

C T5R6PE IS The Tpa\S m I5S!0N COEFFICIENT FC© The may REFLCTING from Th£ 

C 5/6 lNTfrFAC r 4 \ r r c r r L t. C T I NG r R C * ThE <4/3 INTERFACE «ITm 

C PEPRENOICLLAP p;la- I2AT ION- 

75R6PFiTClPEN*T12PEf. •I2 ? REN» T 3«4PEN*TLEP»'N*RS6PEN»TSLPrN»Tu3PFN*(-P 

13‘»PEN»*T*i?P‘'N*TE6PEN*An**ATTC*ATT;»(lATT‘i»AiTS»*«(3.)»*CExP|xS‘4b» 

C T546P* IS ThE T pa *jS“ I S S I ON COEFFICIENT fcp Tur day REFLCTING FROM ThE 

_C - S/6 INTERFACE an: REFElECTINO from. The 4/3 INTERFACE -ITh 

C PARALLEL POLARIZATION 

T c 46PA:TriPA?*T«.PtD.TrTrAR»T34PAD«T4SPAR«R t fcDAR«TELPAR»T47PAR*«-P 

13RPARI*T4SPAC«T'tPAP»ATTl»ATT?*ATT3*«IATTU*ATTSI»Rl3.l»«CEXPIX546» 

XW«44fcl(XEl*XE2^XE3»r.»xE4*irS » / R u 

*4 46:CmPLX(' , .J,-a.«4 6) 

C 

C T 4 4 6 PE IS the TPANSMI5SICN COEFF:cIEN t FC c ThE RAY UJTM ONE INTERNAL 
-t - - REFLECTION IN T h f LlH LAYER bl!) t E R P C NO I CUL A u POLARIZATION 

T446PE: T 31PEN «T 1 CPE f • TZ7PFN*! ?4Pt *;<R4SPE N» < -R J4PENI • Tu 5PEN* T 56PEN* 
lATTl*ATT2*ATT3*CATT4««»i.)>*ATTS*ClXPCX446» 
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IkHbPi IS THE TR A N$M 1 SS T ON COEFFICIENT F0O ThF PAY WITH ONE INTERNAL 

— ®EFlEC TICN IN ThF <* f m IAyEC w 1 Th PARALLEL P0LAOI7ATICN 
T«l4fcPArT11PA 0»Ti;»tC*T;:r'iR*TT«PAl?»ffH^PAV»l-P3HPAP»*T4 5PA(;*T56PIP« 

1ATT1»ATT7*MT *. ( a T T *- •• I 3 . I ) M 71^ «CixP I «4<«6 I 
XW5 56 R I xf l»xE2**E.3«>ifc»3»*XE5 I/PC 

— - XSSfciCKPl X »3. C f - 1-556 > — - 


T 556 PE IS THf TP A NS- I < S I Cf. COEFFICIENT FOPThf »AY WITH ONE INTERNAL 
REFLECTION IN Tnf 5Tn I.AtEP »ITH PERPEND ICULAi? POL API Z A T I CN 

T556PE=T3JPEf. •T12pFN«T25PEf.*Tj4PEN*T«*5P£N»oS6PEN*l-l.,*P*46PEN*T56P 

lEN*ATTl*ATT:*ATTi»ATm*(ATl5*»(3.7»*CE*P|X556l 


T556PA is The TRANSMISSION COEFFICIENT FOP Thf day WITH ONE INTERNAL 
REFLECTION IN ThE 5 Th LAYS 3 WITH PARALLEL POLARIZATION 
T556PA;TC lPAPtT12P*R»T23PaR*T 3‘tPAP*T4SPA»:*R c :6 p AP*|-l. )*P45PAR*T56P 
1AR*AIT1*A1T2*ATT2*ATT4«(ATI5»*I3 . »»*CExP(XS56J 


T255PE IS The TPA\S“ISSICN COEFFICIENT fop ThE pay with Two Internal 

REFLECTIONS In The STh LAYER -ITh PERPENniCUL AO POLARIZATION 

T255PE = T5 56PE «R5 6PEN* C-RN5PEN J»ATT5*A TT5»CE XP»?.*XC5) 

• \ 

1255 PA IS ThE TRANSMISSION COEFFICIENT FOP ThE pay WITH TWO INTERNAL 
REFLECTIONS IN The 5Th LATER -ITh PARALLEL POLARIZATION 
T255PA=T556PE«R5bPAC«(-OHSPAR)*AIT5«ATT5»CExP(2.»xC5) 

• 

T 355 PE IS The TRANSMISSION COEFFICIENT fc® T«r cay with Tw? r E' INTERNAL 
REFLECTIONS IN th- STh layer -ITh RFRPFNDICUL 4= POLARIZATION 
T355P£=T255RE •RScPE*-* ( -RhSPEN »*ATT5* AT T5*CE XP«2.*XC5» 


T 355 PA IS The TRAN EMISSION COEFFICIENT FOP ThE RAY WITH Tncrr INTERNAL 
REFLECTIONS IN the ;Th layer -ITh PARALLEL POLARIZATION 
T355PA;T25jPA*R5t,PAR*t-KH5rAR)*ATTE*ATI5*CEXPI2.*xC5J 


C 

C 


c 

c 


TXSHpF REPRESENT R£ REFLECTION CF RAYS T R J i»pr «T;i4'Pf ♦T®5 t P r F;"“ y/u 

INTERFACE and “STAL E = r u NC PLANE FC® 0 E R PF no I Cl'L A R “ PC L * R I Z A T I ON' 
TX34PrAlTR34PE«TR4:FE-TPEtPE )*R34PFN*T32PEN*T?lPtN»(-l.) » T I 2 P N ♦ T 2 
lJF>£N»MATTl*ATTr*ATT3)**2.l*CEXP(?.*(XCl«xC?*xC3n 

TX34PA S A w E AS T * p E E*C"P T PARALLEL POLARIZATION COMPONENTS 
- . TX34PArtTR34pA*TR4E r A*7?56FA) •03-PAR»T3rPArf*T21PAR*T-l.) *T12PAR*TZ 
1 3R A R • l(ATTl*ATT2*ATT3)»*Z-l*CExP(?.*TXCl«XC?«xC3IJ 

TX45PE S A “E as Txtupi EXCEPT Pf REFLECTION r?C*< n/S T N T E c F A f f 

TX4SPL = eTR3HrZ*TSaSP:»TPS,fef'Ei«RN5PEN* T M3REN*T32 t »EN*T21P£\*C-l.l*T! 

12PCN*TZ3PEN»T3HP£». •((ATT;-ATT2*AIT2*Ani*|»«?.)*cExP(2.*(xCl*XCZ«XC 
2 3 ♦ X C H > I 


C 

_C TX45PA Sa-£ AS T i 3 4 p A EXCEPT F E RE FL E C T I Of. FPj- u/S INTERFACE 

TXN5PA:(TR3<4PA«TR<tSP£*TP56PAI«o««5PaPaTw3PA?«I3‘>PA^«12IFAR«|-I.|*Tl 
12PXR*T23PAR*T3 | *Pah» ( (AIT 1*ATT2*ATTJ*AT74 )»*?.)*CfcXP(2.»( XC1*XC2 # XC 
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2 3*XC4 I I ■ 

. IXSSPf SA-C as TX34PE rxcr?T SERE FLECTION FROM e /6 TnTEpFACC 

TX5lPC:|TPj«F**Tf.4‘rc.IP;6FEl*P'j6 B £N»T54? r *j*T43Pr». • T3?PEN»T?1P« , NM 
1-1. I • T 1 2 p F S • T 3 3 P f *. • T : 4 ; ' ' f. • T 4 5 » «. N • I ( i T T1 • i IT > A T T 2 •ATT4*ATT5l»»2.)» 
2CEXPt2.*fxci*xC2*xC!*xC«*xC5)> 

TX56PA $«“£ AS T * ?4 P A FtCEFT REFLECTION FPC“ S/6 T \Trer»c* 

TX5feP4:tTP34aa«rF4ECA.IO;tp«|. = 5 e) PiP*T64PAP*T4jPA=‘»TJ?PAP*T21PA3*( 

l-l.>»T12?iR»T23PAR*Tj:.PAR*T«5PAR*UATTl*ATT2»ATT3»ATI»i*ATT5»*«2.)* 
2CEXPI2.«lxCl*xC2*xC?*4C«*xCS)» 

tPENOT IS THf TOTiL T B £ » S - I S S I C N C C £ F F I Cl r N T F C THE iUM OF ALL 
RAYS '-ITH P*PPESriCLL43 POLARIZATION 
— lPf KC T: TP-!f. • T C56PF « TC 45 P £ ♦ TR 34PE * T S*. 6 ?E • T4 3 E * T SSfcPE 

1 aTTSSPE *T 3SSPE 

2 «TR 12 PE»T »2 3PE*Tx3 - 4PE»Tif4SPE«TX56PE 

TPARAT IS THE TOTAL TRANSMISSION CCEFFJt I f N T FCR THE SUM CF ALL 
RAYS WITH PAPALufL POLARIZATION 
TPARAT:TPi?*TC56Pi.T345PA»Tc3uPi.T546PAwT4<i6FA.T£56PA 
1 «T2SS ? £ ♦ T 'SSFA 

— . 2*TR12PA*TR2:PA*Tx34rA«TX4SPA«TXS6PA 

TPE REPRE ?fUS' THE “AGNITuCE CF ThE TOTAL PERPENDICULAR T p AN.S*« I SS ION* 
CCtFF ICIEST 

IFERCAbS l TPENOT I . ... 

TPA PEPPESEMS The -ACNITuDE CF THE TOTAL PARALLEL TRANSMISSION 
CCEFF ICIEST 

TPAiCAE' TP AR AT ) . . . 

TC 1° IS The total TPA NS“I r SION COEFFICIENT FOP THE SUM OF ALL 
RAYS «ITH CIRCULAR -01* FIXATION 

TC1RR t T F E N 0 T ♦ TP A R A I 1/2 

TCIRCViCASS ( TCI? ) 


DB PEPRf S r STS A rrci ; . r l F.'C’CR PCP CIRCULAR polapIZA 

ce = :c.<- ii ci t: ; . i 

WRI T£ (e . • I I ’ J , Tt t , ; ; ft , TCI"C V , 0? 

50 FORMAT ( S-a,F 1C ',5x,F5.3,SX,FS.3,SX,F5.3,5X f r 7.3l 
20 CONTINUE 

IFU.EO.il GO TO 3 
END 


ppl*PIZAT!ON 




